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A future measurement of cosmic microwave /i-distortions by an experiment with the specifications 
of PIXIE will provide an equivalent 3a detection of the running of running of the spectral index of 
scalar perturbations, 13 = da/ din k, ii fi > 7.75 x 10~*, covering much of the PIXIE sensitivity range. 
This corresponds to a resolution limit of /3 > 0.015 which is relatively large given any presumption 
of slow roll, a result of the current tight constraints on a < on CMB scales. We show that a 
detection of /3 at this level is in conflict with slow roll conditions if the primordial signal can be 
distinguished from any post-inflationary contamination. 
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I. INTRODUCTION 

Understanding the character of the primordial density 
perturbation has been a chief endeavor of modern preci- 
sion cosmology. Our knovifledge of the shape and statis- 
tics of the primordial power spectrum has grown substan- 
tially over the past two decades, from the Cosmic Back- 
ground Explorer's first glimpse of the cosmic microwave 
background (CMB) anisotropics, to today's swarm of 
ground-based and balloon-borne CMB and large scale 
structure observatories. Our knowledge of the primor- 
dial density perturbations derives almost exclusively from 
these data sources - the temperature and polarization 
anisotropies of the cosmic microwave background and 
galaxy correlation measurements - both of which probe 
large length scales. These observations have revealed a 
nearly scale invariant spectrum of Gaussian, adiabatic 
perturbations on comoving scales k < OA Mpc^^, in 
agreement with the simplest models of inflation. 

But inflation might not have been so simple. Devia- 
tions from Gaussianity and adiabaticity are possible, and 
the power spectrum might deviate from a power law |TJ- 
[5] . A spectrum with a nontrivial scale dependence might 
indicate a complicated inflaton potential exhibiting non- 
slow roll behavior. A correspondence exists between the 
scale dependence of the power spectrum, P(k), and the 
fleld-dependence of the inflaton potential, V{4)): specifi- 
cally, the n'^'-order term of the Taylor expansion of V{4)) 
in (j) is linear in the n'^'-order term of the Taylor expan- 
sion of P{k) in Infc 4J. By resolving the scale depen- 
dence of P{k), in particular the higher-order runnings of 
the spectral index: a = dns/dlnk, j3 = da/dlnk, etc., the 
paradigm of single field, slow roll inflation can be tested. 

In order to observe these higher-order runnings, the 
spectrum must be measured across a wide range of length 
scales. The Wilkinson Microwave Anisotropy Probe 
(WMAP) and Sloan's observations of baryon acoustic os- 
cillations combine to probe scales 10""* Mpc~^ ^ ^ ^ 0-1 
Mpc~^; while these data exhibit a preference for a < 0, 
it is not statistically significant [5j. The Lyman alpha 
forest, which probes large scale structure on scales /c « 1 



Mpc~^ [Ms], is promising, but converting flux measure- 
ments to a density power spectrum is difficult with our 
currently incomplete understanding of the properties of 
the intergalactic medium. Recent ground-based observa- 
tories, like the South Pole Telescope[^nd the Atacama 
Cosmology Telescope]^ have begun to place nearly de- 
cisive limits on the running of the spectral index by 
probing the damping tail of CMB fiuctuations: a — 
-0.034 ± 0.018 and a = -0.024 ± 0.0112 |10j, and 
Planck will do even better. Current data, however, can 
go no further: in order to resolve higher-order runnings, 
we need a longer lever arm. By extending the observa- 
tional window to smaller scales, P{k) can be recovered 
across a greater range of In/c, and subtle higher-order 
terms in its Taylor expansion will become relevant. One 
promising probe of small scale fiuctuations is the 21cm 
transition of neutral hydrogen, which will extend the ob- 
servational window down to scales k ~ 100 Mpc~^ |11| . 
and resolve a at the level of 10"'* [HIS]. We can expect 
these results in the coming years. 

But it might be possible to go even further: the ther- 
mal spectrum of the cosmic microwave background radi- 
ation traces fluctuations on scales 50 Mpc~^ k < 10^ 
Mpc~^, offering an unprecedented glimpse into the char- 
acter of the small-scale power spectrum. This is a fun- 
damentally different kind of measurement: rather than 
probing the fluctuations themselves, the focus is on the 
dissipation of these fluctuations into the baryon-photon 
plasma and the resulting deviation of the thermal spec- 
trum of the radiation field from that of a black body. The 
proposed PIXIE [T3] mission will be sensitive to these 
spectral distortions at a level useful for power spectrum 
reconstruction, and there has been recent interest in its 
prospects for improving constraints on a |15H18j . In this 
paper we investigate the potential of such a measurement 
for resolving the running of the running of the spectral 
index, P = da /dink. 



^ http://polo.uchicago.edu/ 

* Email ^brian.powell007@gmail.com| ^ http://www.physics.princeton.edu/act/ 
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II. CMB DISTORTIONS 

Fluctuations in the baryon-photon plasma below the 
Jeans scale oscillate as acoustic waves during radiation 
domination. The energy of these fluctuations is dissi- 
pated via diffusion damping, resulting in a redistribu- 
tion of this energy from the fluctuations into the baryon- 
photon plasma [TUH2S]- At early times (z > 2 x 10®), 
efficient double Compton emission and bremsstrahlung, 
and the up-scattering of these photons by energetic elec- 
trons quickly thermalizes the baryon-photon plasma, pre- 
serving its black body spectrum. As the universe cools 
- between redshifts 5 x 10^ ^ z < 2 x 10® - Comp- 
ton scattering ensures statistical equilibrium of the radi- 
ation field, but the production of photons is insufficient 
to maintain a black body: a frequency-dependent chem- 
ical potential, /x(aj), develops and the spectrum shifts to 
that of a Bose-Einstein distribution. So, while small- 
scale perturbations are erased before recombination, they 
leave behind a record of their existence: the energy lost 
by the perturbations is injected into the CMB distort- 
ing the spectrum away from a black body by an amount 
that depends on the shape and amplitude of the power 
spectrum on the relevant scales. 

The time evolution of the /z-distortion amplitude obeys 
[251 
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where tjjc is the double Compton thermalization time, 
and Q / p~f is the fractional rate of energy dissipation, and 
double Compton emission has been assumed the domi- 
nant source of photon production [551 HZ] ■ This has the 
solution 
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where = 5 x 10"* and zdc — 1.98 x 10®. The energy 
density of the acoustic waves is Q = 3p-yC^((5^(x)^)/4, 
where w 1/3 is the squared sound speed, and the pho- 
ton density perturbation, 
age. 



(5^(x), has the ensemble aver- 



(<5^(x)2) = 



d^k 

W) 



(3) 



where P{k) = A^P'(fc), the transfer function, and 
P^(k) is the primordial spectrum generated by inflation. 
Since we are considering acoustic waves below the Jeans 
scale during radiation domination, the transfer function 
is appropriate to modes that are well inside the horizon. 
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where is the sound horizon and the diffusion scale, fc^i, 
varies with redshift during radiation domination as 
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where = 8 (l35iJo\/n7neocrT) = 5.92 x 10 
Mpc"""^. Here rieo is the free electron number density 
at z = and ar is the Thompson scattering cross sec- 
tion. The portion of the power spectrum, P'(k), that 
is primarily responsible for generating the /^-distortion is 
set by the diffusion scales at the beginning and end of 
the epoch, 5 x 10^ ^ z < 2 x 10®. Equation ^ gives 
k e [50, Mpc-i. 

With these results, the amount of //-distortion Eq. ^ 
can be well-approximated by |18| 
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where cos^(fcrs) has been replaced by 1/2 - its average 
value over one oscillation. In contrast to the anisotropy 
measurements for which CMB data on scale £ constrains 
the spectrum at corresponding wavenumber, k, the fi- 
distortion measurement comes in the form of a weighted 
integral constraint on P*(fc) over the /x-distortion epoch. 



III. OBSERVING (3 

The question of whether an observation detects a phys- 
ical quantity, 0, is a problem in hypothesis testing: are 
the data sufficient to overturn the null hypothesis, in 
which 6 = 07 In Bayesian statistics, hypotheses are in- 
terpreted as models, and hypothesis testing proceeds via 
a comparison of the Bayesian evidence of the first model, 
p{d\A4o), with that of the second, p{d\Aii), 



Bn 



pid\Mi) 



(7) 



where d is the data and Ai the model. The Bayesian 
evidence is the average of the likelihood function, 
C{8\d, M), weighted by the prior probability of the pa- 
rameters, 9, 



p{d\M)= C{0\d,M)Ti{e\M)de. 



(8) 



Once the evidence of each model is computed, the Bayes 
factor Eq. ^ can be evaluated and inferences made 
based on the Jeffrey's scale [55]: InBoi < 1 means that 
the data is insufficient to distinguish the models, while 
1 < lni?oi < 2.5, 2.5 < lni?oi < 5, and IuBqi > 5 signify 
positive, moderate, and strong evidence, respectively, in 
favor of A^o- In what follows, we will quote results that 
satisfy the threshold for strong evidence in favor of A^i, 
I In Soil > 5, corresponding to a 0.993 posterior proba- 
bility for the favored model 29J. 

We are interested in the prospects of detecting and 
constraining higher-order runnings of the scalar power 
spectrum, specifically, the running of the running of the 
spectral index, /3 = da/ dink, with future observations 
with the specifications of the PIXIE experiment ^3], a 
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FIG. 1. Marginalized posteriors of a for different fiducial 
amplitudes of /^-distortion for model Mo- The value of the 
best fit Xe// = — 21n/I is shown for each case. 



proposed nulling polarimeter for the measurement of pri- 
mordial gravitational waves. The model that we wish to 
study, A^i, has a scalar power spectrum given by 



(9) 



The base model, A^O: has /3 = and so is contained 
in A^i as a special case. This is a standard exercise in 
Bayesian model comparison: is the data sufficient to war- 
rant the addition of the parameter /3? Equivalently, the 
model comparison determines whether the likelihood of 
the model fit to the data is sufficiently improved, beyond 
statistical coincidence, when the new parameter is in- 
cluded. Because the models M.o and M.i are nested, we 
make use of a representation of the Bayes factor known 
as the Savage-Dickey density ratio |30j . 
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where p(/3|d, A^i) is the marginalized posterior distribu- 
tion of /? and 7r(/3|A^i) its prior probability under model 
Ml, and both are evaluated at /? = in Eq. ([TOl. 



The Savage-Dickey density ratio is a convenient route 
to the Bayes' factor that does not require an evaluation 
of the Bayesian evidence, which typically involves a high- 
dimensional integration that can be computationally de- 
manding. It is nonetheless exact for nested models with 
separable priors. 

Current data is consistent with a null detection of /3, 
and so our objective is to determine how large the /i- 
distortion must be in order to positively detect /3 with 
PIXIE, i.e. to substantiate its inclusion as a parame- 
ter in the model. We make use of two data sources: real 
cosmic microwave background temperature and polariza- 
tion anisotropy data, and a notional /^-distortion signal 



of the quality expected from PIXIE. The former data 
set includes the latest measurements from the Wilkinson 
Microwave Anistropy Probe (WMAP) [S] and the South 
Pole Telescope (SPT) [HT. These data are well-described 
by the standard inflationary cosmology with base param- 
eters: the baryon and cold dark matter densities, Q,i,h? 
and flch^, the ratio of the sound horizon to the angular 
diameter distance at decoupling, 9s , the optical depth to 
reionization, r, and a power law spectrum of density per- 
turbations parameterized by amplitude. As , and spectral 
tilt. Us- The SPT data requires three additional fore- 
ground parameters: the Poisson point source power from 
randomly distributed galaxies, D^^, the clustered point 
source power, Df^, and the secondary emission from the 
Sunyaev-Zeldovich effect from clusters, Aa,z- The model 
likelihood calculated from these parameters, which we 
denote Cce to indicate that it refers to the WMAP and 
SPT Ci anisotropy spectra, is described in O [10] . 

The second data source, that arising from fi- 
distortions, is notional. For a given fiducial amount of 
/i-distortion, /ip, the likelihood of the model with a theo- 
retical value of 11 given by Eq. ^ is computed assuming 
a Gaussian posterior for fi, 



Cf, cxp(/i)/7r(/i). 



cx exp 
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where the prior, 7r(/i), is taken to be flat within the range 
fi e [0, 10"''], and ct^ — 10^^ Then, for a particular 
set of parameter values, the full likelihood of the model 
can be determined by simply summing the log-likelihoods 
of the two data sources. 
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To summarize, we are considering an overall likeli- 
hood that derives from two separate data sources, CMB 
anisotropics and /i-distortions, that probe different cos- 
mological scales. 

Since the real CMB anisotropy data are fixed, our task 
is to vary ^.q across the prior range until we achieve 
I In Boil > 5. This process sets up an interesting ten- 
sion between the two data sources. On the one hand, 
we have CMB temperature and polarization anisotropy 
data that is well-fit by the model TWq with /? = 0. In 
fact, the simpler model with a = is still marginally 
preferred by current data [32j , and is perfectly consistent 
with an amount of /^-distortion below the PIXIE detec- 
tion threshold oi ^ ^ 10~^. However, as the fiducial 
amount of /i-distortion is increased, even the model with 
running, Mq, struggles to provide an acceptable fit. This 
is illustrated in Figure 1: as /iq is increased, a must get 
larger to accommodate the increase in small scale power 
needed to generate the fiducial amount of /^-distortion. 
But Mo can only provide increased small scale power at 
the cost of simultaneously increasing large scale power, 
and in an effort to fit ever larger values of /iq, the model 
is pulled out of agreement with large scale CMB data. 
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This effect signals the need for an additional parameter, 
P = da /dink, that gives a spectrum with increased small 
scale power - to accommodate the large /^-distortion - 
without strongly modifying the shape of the large scale 
power spectrum that is well-constrained by current CMB 
data. 

In order to calculate Bqi for a particular value of /io, 
the posterior, p{(3,d\Mi), and prior, 7r(/3,A^i), distri- 
butions must be determined. The prior encodes our 
current knowledge of /?, and is taken equivalent to the 
posterior distribution of /3 obtained from WMAP-I-SPT 
data: it is close to Gaussian with mean /3 — 0.017 
and standard deviation cr^ = 0.024. The poste- 
rior distribution is obtained using Markov chain Monte 
Carlo (MCMC) over the parameter space of model Aii: 
{Qbh^ n,h^, es,T, A„ n„ a, /3, r, D^, Df^,Asz}- For 
each sample, Eq. ^ is numerically integrated with P'^{k) 
given by Eq. ^ to obtain the //-distortion and the over- 
all likelihood, £tot, determined from Eq. (12). This pro- 
cess is iterative: if the i?oi thus calculated for the chosen 
fiducial value of /to does not satisfy | ln_Boi| ~ 5, a new 
fiducial /xq is selected and another MCMC run is begun. 



In Table 1 we present the results of our analysis for 
different values of /to in the range jio G [10~^, 10"'']. For 
each value of /to, six chains were run using the parameter- 
ization of model A^i, and chain convergence was consid- 
ered achieved when the Gelman-Rubin statistic satisfied 
R~ 1 < 0.05. The CosmoMCj^ software package [33] was 
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TABLE I. Bayes factor, | In Soil, indicating varying levels of 
evidence for /3 7^ for different fiducial /t-distortion ampli- 
tudes, /IQ. 

used for sampling and chain analysis, with modified like- 
lihood routines to calculate £tot- We find that PIXIE 



http:/ / cosmologist.info/cosmomc. 
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will detect /? 7^ if /Lto > 7.75 x 10~*, covering much of 
the PIXIE sensitivity range (/z > 10~®). The correspond- 
ing resolution limit of the running of running is relatively 
large: /3 > 0.015. This is principally a result of the good 
constraints on a < on CMB scales: /3 is forced to be 
big to yield a large enough /^-distortion to be detected by 
an instrument like PIXIE. 

In addition to detecting /3 7^ 0, PIXIE will place rel- 
atively good constraints on its value: (3 = 0.015io gi 
at 2a, for fiQ = 7.75 x 10^*. Current constraints on 
the spectral parameters for WMAP-I-SPT are compared 
with projected constraints when PIXIE is included at 
this level in Figure 2. It is interesting to note that while 
/i-distortion data greatly improves constaints on j3, there 
is little change in the posterior distributions of either 
or a compared to current results. This is because both 
Us and a are already well-constrained by WMAP+SPT: 
the existence of sufficient small scale power to generate 
a detectable level of /i-distortion establishes the need for 
/3 and constrains it without affecting the tight limits on 
n» and a on CMB scales. 



IV. IMPLICATIONS FOR SLOW ROLL 
INFLATION 

The minimum amount of running of running that can 
be detected with statistical significance (strong evidence, 
or 0.993 posterior odds) is relatively large: /3 > 0.015, 
possibly in conflict with slow roll inflation, provided that 
the observed /i-distortion is primordial. The slow roll ap- 
proximation holds whenever the acceleration of the scalar 
fleld, 0, can be neglected relative to the cosmological 
drag, (j> ^ 3iJ0, with the result that 



3H(j) « -V', 



(13) 



where dots denote time derivatives and primes denote 
derivatives with respect to 0. This is the case as long 
as the scalar fleld is potential-dominated: 0^/2 <C V {(/)). 
Slow roll can be maintained as long as the potential is 
sufficiently flat; this condition is often stated in terms of 
the smallness of the parameters. 
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where the dots denote terms involving higher derivatives 
of H{(j)). But how many terms must we consider? How 
small do they need to be? The slow roll approximation 
can be formally derived by considering the Taylor expan- 
sion of the Hubble parameter [31] , 
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where Bi = \J eo/2 and the subscript "0" indicates that 
the coefficients should be evaluated when the pivot scale, 
fco, exits the horizon, and we have chosen = 0. The 
coefficients, 
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generate the hierarchy of slow roll parameters derived 
from H{(j)): ^Xh = ??, ^Xh — S,^, etc. In this representa- 



tion, a cosmological solution given by Eq. ( 14 ) truncated 



at some order M is interpreted as a perturbation expan- 
sion about de Sitter space. 

The parameters e and the ^Xh determine approximate 
expressions for the spectral parameters: rij, a, /3, 
that are themselves coefficients of the Taylor expansion 
of InP(fc) in Ink. The lowest order spectral parameter, 
ris, is lowest order in slow roll, ns(e, r/, • • • ), where ••• 
includes all higher-order terms e^,er],r]'^, and so on. The 
second order spectral parameter, a, is likewise second 
order in slow roll, ck(^^, e^, ery, • • • ). The slow roll approx- 
imation can be considered valid if the n*'*-order spectral 
parameter is dominated by the n*''-order slow roll param- 
eters, with higher-order terms comprising ever smaller 
corrections. It is therefore a prediction of slow roll in- 
flation that higher-order spectral parameters should be 
vanishingly small, and that a general solution should be 
described by only the few lowest-order terms. 

We begin by writing the spectral observables as follows 



= 1 - 4e + 2// - 2(7^^ _ ^(^(j _^ ^ _^ lOC)e?7 
a = -2^2 _ -^Qg^ _ _^ UC)ef 

+2Cr]f + 2C'^Xh 
(3 = -Ue^^ + 2tj£,^ - 2C£,^ + 2 - 2C ^Xh 
-32e^ + 626^7/ - 20er/2 _ (46 + hQC)e^i^ 
+ {2& + 48C)e77C^ - 'iCifS,^ + (10 -I- 20C)e ^A^ 
-6C?7 ^A^, (16) 

where C ~ —0.73. The subscript "0" has been omit- 
ted for brevity; it should be understood that the above 
expressions are to be evaluated at the pivot scale, /cqj 
and associated field value, 0o- We have written the ob- 
servables to next-to-leading order in the slow roll param- 
eters because we would like to study the contributions 
of these higher-order terms, especially for /3. We first 
set ^A// = 0, truncating the hierarchy at order M = 3. 
We then seek a cosmological solution, parameterized 
by {e(</)), r?(0), ^^(0), "^A/f ((/>)} that yields observables in 
agreement with current data in combination with a dis- 
covery of /i-distortion by PIXIE with /io = 7.75 x 10~^ 
(cf. Figure 2.) 

In addition to satisfying parameter constraints, the so- 
lution must yield sufficient inflation to solve the horizon 
and flatness problems: A^ > 30 e-folds if inflation is to 
occur above the electroweak scale, and A^ w 55 (60) for 
GUT (SUSY) scale inflation. The number of e-folds of 
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FIG. 3. (a) Distribution of models generated stochastically at third order in slow roll. The models are in agreement with current 
CMB data, but they fail to provide a sufficient number of e-folds, A*', of inflation or yield large enough /3. (b) Distribution of 
models generated stochastically at fourth order in slow roll. While most models still provide insufficient inflation, many are 
successful with large /3. 



expansion is given by 
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where H{cj)) is determined by the slow roll parameters 
up to order M = 3 via Eq. ( Il4|. For a given set of 
parameters {eo,'7o,Co' ^^ho} Eg. (|17[) can be integrated 
numerically to obtain the number of e-folds of expansion 
that occur after the scale fcp leaves the horizorjj We 
stochastically sampled the e — 77 — — ^Xh space, and 
integrated Eq. (17) for each set of initial conditions. We 



collected solutions with r < 0.3, e [0.94, 1.01], and a G 
[—0.08,0], ranges that approximate the 95% confidence 
limits obtained on the parameters in the previous section. 
No constraints were placed on TV or /3 so that the full 
space could be explored. As we show in Figure 3 (a), 
not only do the models fail to provide enough inflatiorj^ 
they all give (3 < 0.005 - just outside the 2a limit. This 
result shows that if /3 is controlled only by leading-order 
terms in the slow roll expansion, it is not possible to find a 
working single field inflation model that flts observations. 

What about the next-to-leading order term, ^Ah? 
If we truncate the slow roll hierarchy at order 
M = 4, our analysis shows that not only is 
sufficient infiation obtained, but /3 can now be 



* This is equivalent to solving the system of flow equations trun- 
cated at order M = 3 [381139] , 

^ A similar investigation was carried out in 1401 regarding the im- 
plications for slow roll inflation of a near the WMAP3 centroid, 
a ^ -0.05. 



large (Figure 3 (b)). One such model is given 
by the parameter values: {eoi^iCo' ^•^Hi ^^h} ~ 
{0.015, -0.01, 0.035, 0.018, 0.0025}. Not only are the first 
four parameters of equal magnitude, but the fifth param- 
eter, '^Xh, while only amounting to a small correction ac- 
cording to the slow roll approximation, strongly affects 
the inflationary dynamics and allows for ^Xh to become 
large. Our analysis therefore shows that for /3 > 0.005, 
the slow roll expansion must be taken to at least fourth 
order to obtain viable solutions, many of which are de- 
scribed by slow roll parameters of roughly equal magni- 
tude and relevance. Given this observation, there is little 
reason to expect that higher-order terms can be safely 
neglected. We conclude that a statistically significant 
detection of running of running by PIXIE, (3 > 0.015, is 
in conflict with single field, slow roll infiation. 



V. CONCLUSIONS 

The discovery of higher-order runnings of the scalar 
power spectrum has important implications for the na- 
ture of the inflationary era. Large higher-order runnings 
threaten the validity of the slow roll approximation, and 
disfavor an explanation in terms of single field, slow roll 
inflation. Current data from WMAP-^SPT have almost 
conclusively detected a ^ 0, and Planck will furnish a 
detection if the value of a is close to current best-fit es- 
timate, a w —0.024. Meanwhile, the running of run- 
ning, l3 — dhia/dlnk, is not well constrained because 
today's observatories do not probe the primordial spec- 
trum, P{k), across a sufficiently large range of length 
scales to resolve it. In this work, we investigated whether 
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a future measurement of CMB /i-distortions by an ex- 
periment like PIXIE, taken in combination with CMB 
anisotropy data, would provide a sufficiently long lever 
arm to resolve these higher-order runnings. 

We framed the detectability of /3 as a problem in 
Bayesian model comparison between two power spectra: 
one with j3 and one without. We computed the Bayes 
factor using MCMC for different fiducial values of /io in 
order to determine how large the notional /x- distort ion 
signal needs to be to favor the running of running model. 
We find that if /xg > 7.75 x 10~®, the inclusion of 13 in 
the model is required by the data with strong significance 
(equivalent to 0.993 posterior odds.), covering much of 
the PIXIE sensitivity range. The best-fit value associ- 
ated with detection is /3 w 0.015, which is relatively large 
given any presumption of slow roll. This is a result of 
the tight constraints on a on CMB scales: with a < 0, 
j3 must be moderately large and positive in order for the 
spectrum to produce enough /i-distortion on small scales 
to be detectable by PIXIE. Constraints on /3 are also 
improved by around a factor of 5 over current estimates. 

We next studied the implications that large (3 has 
for single field slow roll inflation under the assump- 
tion that the primordial distortion is distinct from post- 
inflationary sources. We performed a flow analysis, fo- 
cusing particularly on the next-to-leading order terms in 
the slow roll expansion and their role in the inflation- 
ary dynamics and observables. It was found that models 
taken to leading order [M — 3) failed to yield sufficient 
inflation or satisfy constraints on /3, while a truncation 



at next-to-leading order (M = 4) gave models with very 
different inflationary trajectories that succeeded in meet- 
ing these criteria. The inclusion of ^\h, while amount- 
ing to only a small correction according to the slow roll 
approximation, strongly alters the inflationary dynam- 
ics. Successful M = 4 models had slow roll parameters 
of roughly equal magnitude and importance, challenging 
the slow roll prescription that higher-order terms in the 
expansion should indicate ever smaller corrections to the 
inflationary observables and trajectories. We conclude 
that a detection of running of running by PIXIE will be 
in conflict with single field, slow roll inflation. 

The search for CMB /x-distortions is a promising av- 
enue for challenging the prevailing conception of infla- 
tion. It is important that work continues on understand- 
ing the possible foregrounds that might contaminate any 
^-distortion signal, as well as identifying possible post- 
inflationary sources of /x-distortions that might mimic 
a primordial signal. If a measurement of primordial /x- 
distortions by PIXIE results in a detection of /3, then this 
will be a strong indication that inflation either exhibits 
non-slow roll dynamics or is driven by multiple fields. 
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